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™> CT »TfSniI ^PPn« n; T.T.TTTn STOB-^E SYSTEMS 
r»w» BB PAfprpnce Related-"Applicat ions, 
This is a continuati^i^-part of serial no. 
09/083,732, filed Mav^l998. which is a continuation-in- 
part of applicatj^erial no. oa/SSa, 129, filed, June. A, 
1996 now U^Patent No. 5,757, ff«T4, and a contrinuation-xn- 
part'o^Trnational patent application no. PCT/TOST/OWW. 
^l_ ^duuLi un:i l f ili n g- dam u £ J uaeJ^y 

RArVarQUpH of the TnvBntion .. 

The invention relates to monitoring the operational 
performance of fluid storage systems . ^ > 

Large quantities of liquids and similar materials 
are often stored in bulk, storage containers .or tanks, which 
may be located above-ground, partially above -ground , or 
completely below ground. Such, cor^inexs or tanks ar*. 
generally connected by piping to flow-meters or dxspensers. 

For example, underground, storage tanks (UST's) and, 
occasionally, above-ground storage tanks (AST's) are used to 
store petroleum products and fuel, to be dispensed at 
automobile service stations, trucking terminals, automobxle 
rental outlets, and similar operations through gasolxne, 
diesel, or kerosene dispensing, pumps . Fuel product xs 
generally delivered to such facilities by a gravity drop 
from a compartment in a wheeled transport means such as a 
fuel delivery truck or an introduction of product through an 
underground piping system. AST's or UST's are often located 
at central distribution locations so that product can be 
subsequently withdrawn from the tank system to be 
transported for delivery to a variety of such facxlxtxes. 
distribution location with UST's or AST's may receive 



10 



15 



liveries of product f—. e.g.. a pipeline spur, wfcaeled 

k and storage containers- is difficult or. 

such tanks ana stor^ amount 
• «oo-iMe The various methods for identifying trie 

repeatability, and performance ^^p™^ 

^vices which measure « ^ f ^f^lr. — « 
ch e storage containers ^ dx f g- ^ 
mav not be temperature compensated. xne 

in Slivered to the tank system is often measured 
actually delivered to t« Rather, the owner 

inaccurately and, frequently, not at all. Rather, 
To^rato/of the tank or vessel ^^^^ 
voiced amount of product ^"^^^^^ of 
introduced to the tank system, without having y ^ 
confirming whether the invoiced amount of produ 
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25 



" C ° rr cogently, effect^ ^ 
,^mr»i ieated by the numerous errors in trie 

laseline for management, planning and decisionmaking. 

Kffective management *^*J^^°™i *» 

1 Accurate measurement or trie voxu 

system. ^ . nation of ttl e volume dispensed from 

2 Accurate determination octxi 

predetermination of the amount of product 

" fhe tank system — — «~ 

s . Rapid identification of leakage from the tank 
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system . 
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s con t iruaou S :mani.t,oring and. diagnosis of the operating, 

performance of all of the cotnponent measuring 
devices of the system. 

continuous., anaiysis at sale* dac*- .to predict demand 
5 of product from the system. 

8 Determination of optimal reorder times and 

quantities as a function of ordering, 
transportation holding, and penalty costs in order , 
to minimize total costs : of operation, and/air t» 
L0 maximize profits. 

Traditionally, these functions were performed 

«^t- »t- all Volume measurements 
crudely or, in many cases, not at all. voj. 

were, and in many instances still are. based -."T^f 
Pledge of the geometry, dimensions, and <~**«»^£ 
15 the storage vessel. Also, dispensing meters are frequently 
miscalibrated. This is true even when tan* systems^ are 
reflated, due to the breadth of tolerance perked, for 
individual sales as related to total tank volume. For 
example, deliveries from the delivery vehicle are almost 
20 always unmetered. additions of product from defueling 
vehicles are typically undocumented, and theft of the 

oroduct is not uncommon. 

Leakage of product has, in recent years, assumed a 
dimension far in excess of the mere loss of the P-^«. 
2 S Environmental damage can. and freouently does, expose the 
operator to very large liabilities from third party 
litigation in addition to U.S. Environmental Protection 
Agency (EPA) -mandated remediation which can cost in the 
range of hundreds of thousands of dollars. The EPA s 
30 recrements for leaK detection are set forth in EPA 

• „u*- Tallr on Tnn V-°- T,e>ak De fection 
TJn 510-K-95-003 , <3 1- -r-a i ahf TalK OH ianr.=> 

(July 1991) « ^rir W^W^"" r n f.rprv-y 
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To address these concerns,- Statistical Inventory - 
•lT at ^T(SIR) was developed. The SIR method consxsts 
Reconcxlxatxon (SIR) was *- HfiM a n of the 

of a computer-based procedure "* ch of cto 

sou- of error noted ~ ^-^^HLo the 
various and unique patterns that are at cumulative 
inventory data and, in partxcular, into th 
variances in the data when viewed as functions P 
height, sales volumes, and time. - 

« ,,-T-y of rhp Tnventxori, . 

. «^ n > n t-n an automatic SIR. 
The present invent! on/ relates to an au 

Th e invention addresses a variety of ^^t^. 

operational and environmental issues associated 

b ul* storage of lipids or pourable ^ ^ 

The present invention is an application o 
greatly enhances the ability to manage a_ 
effectively. It provides the means to *aracterr 
the geometry, dimensions, and =°<^"° delive ries and 
vessel, identify overages and shortage s in ^ 
unexplained additions and removals of P rod "«' P 
an accurate assessment of overall dispensing meter 

calibration. In addition, by £ ~ ^ tif tca tion 

*-v.-. nreaent invention permits iaeiaiJ - J - 
rXr^-C* -I ^llon per hour in all of 
of leakage at rate* . _ e By increasxng 

its estimates to any P~-ibed ~ e » • ^ ^ 
the number of measurements taken, 
derived at any desired level of tolerance. 

The method of :the present invention makes no 
The metnoa measurxng 
assumptions as to the precxsxon of ^ of w . ion 

^ in various system conf xguratxons . 
devices used m ^ £rom the data alone, 

and, calibration accuracies are u=i 



Also it is not assumed that the Xanic. system is leak free; 
the leak status of the system, is determined from the aata 

alone. . , 

The method derives tank geometry, dimensions , and . 

configuration, and their impact on the totality of 
cumulative inventory variances, as a function of product 
height in the tank. Correctness of dispensing meter 
calibration is verified in a similar manner by testing for 
randomness of cumulative variances- as . a function of varyxng 
sales volumes. Having confirmed that such- remaining 
residual variances are random, -reflecting only the inherent 
random noise of the measurement devices, the present method 
analyzes departures of the cumulative variance from the 
bounds determined by the calculated random noise level. All 
calculations as to the volumes added, removed, metered, or 
leaking are based upon extended successive, simultaneous 
observations of meter and gauge, readings. The nuxnber.of 
observations incorporated in each such calculation is 
determined by computing confidence bands for. the parameters 
of interest and extending data collection as necessary to 
achieve predetermined tolerances. 

For example, the method of the present invention is 
capable of distinguishing between continuous losses 
consistent with leakage and one-time unexplained removals of 
the fluid product from the tank. The method may be used to 
ensure the accuracy of computed, delivery volumes, which are 
determined and reported with confidence boundaries 
calculated for estimated delivered quantities. 

The method can also be used to control and monitor 
the accuracy of purchase- costs of fluids such as petroleum 
which are delivered to. tanks. For example, motor fuel 
retailers may be charged by wholesalers for either net or 
gross. volumes purported, to have, been delivered. * 
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15 



20 



25 



determination that purchase.charses are a PP 
requires frequent simultaneous readings =£ sales_ tank 
volumes and temperatures, which can be accomplished using 
fhe method of the present . inventxon. • 

to accomplish these goals, the present invention 
involves estimating changes in product volume in a tank 
D ased on multiple data points and their T 
errors measured continuously over , a peridot . t*—. ^ 
software program is used to implement an "-^^f» 6 
employs concepts from matrix theory and »^ 

s ra t istics. The ^^^-^^z^rr 

of a matrix and its transpose by successive a 
partial products of partitions of the matrix and their 
partial f „,._w o-rtitions to minimize the 

corresponding transposed matrix P""" compE essed 
storage requirements of the data collected The 
macrii data constitutes a complete and sufficient statistic 
matrix; oata con The algorithm thus- permits 

for the parameters of interest. The aigo 
he ^cumulation and storage of a. large amoun t ~ * 

condensed form without sacrificing ^™?£^Lt 
information, to obtain a statistically signet result 

with the required accuracy a| ^^^^ft^s a 
In general in ^ = aspect ^ / system 

m ethod of monitoring a fluid storage an ^sp ^ 
including a dispensing apparatus . A quantity o 
alsPensed from L system using the dispensing «P-»™ 
eased on an authorization. A plurality of -surement data 
ls collected from the ^^T"^^ — 
C^redT a^ P-lity of measurement 

xs store ^ d to calculate a volume of 

data, are statisticaxj-y auo-uj collected 
fluid based on the plurality of measurement data collected 

from the dispensing apparatus. 
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implementations of -the •.imrentioa may also include 
one or more of the following features. The authorization 
may be based on information relating to a user requesting 
that the quantity of fluid be dispensed, from the system. 
5 The authorization may be provided by a fuel access control 
system. 

The method may further include transmitting the 
plurality of measurement data, to a processor to perform the 
statistically analyzing step. The • plurality of . measuxemerrt 

10 data may be stored in a compressed matrix format, whxch may 
be a product of a data matrix and a transpose of the data 
matrix. The dispensing apparatus may include a totalizer 
and a meter. *.***n//*f c&'^J, 

in general, in another aspect, the invention 

15 features a method of monitoring a fluid, storage and 

dispensing system including a plurality of dispensing 
apparatus. A quantity of fluid is dispensed froa the system 
using the plurality of dispensing apparatus, baaed, on an 
authorization. Measurement data is simultaneously collected 

20 from the plurality of dispensing apparatus in a form 

readable by a computer. The collecting step is repeated to 
obtain a plurality of the measurement data, which are. stored 
in a memory. The stored plurality of measurement data xs 
statistically analyzed to calculate a volume of fluid based 

25 on the plurality of measurement data collected from the 
plurality of dispensing apparatus. 

implementations of the invention may also include 
one or more of the following features. The authorization 
may be based on information relating to a user requesting 

30 that the quantity of fluid be dispensed from the system. 

The authorization may be provided by a fuel access control 
system. 
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The method may further, include, transmitting the 
plurality of measurement data to a processor to perform the 
statistically analyzing step. The plurality of measurement 
data may be stored in a compressed matrix format, which may 
be a product of a data matrix and a transpose of the d * ta ^ 
matrix. ^ .cAt,**^ 

In general, in another aspect, the invention £ 
features- a fluid storage and dispensing system including a . 
tank for storing fluid and. -a dispenser far dispensing fluid . 
from the tank. A fuel access control system enables the 
dispenser based on an authorization. Measurement apparatus 
collects measurement data corresponding to a quantity of 
fluid dispensed from the tank by the dispenser. A processor 
statistically analyzes the measurement data to calculate a 
volume of fluid in the system. 

Implementations of the invention may also include 
one or more of the following features- The system may 
further include a memory for storing the collected 
measurement data. The measurement apparatus may include a 

totalizer and a meter. 

The fuel access control system may be activated by a 
coded medium which contains identification information 
pertaining to the user and which constitutes a request to 
dispense fluid from the system. The authorization may be 
based on information about a user requesting to dispense 

fluid from the system. 

The processor may store the measurement data in a 
compressed data matrix format, which may be a product of a 
data matrix and a transpose of the data matrix. 

Ah advantage of the present invention that the 
accuracy and consistency of devices used to measure volume 
of product added to, removed from, and present in a fluid 
storage system may be determined. 

- 8 - 



Another advantage of .ther present , invention .is that . ■ 
additions of material to the system, but not recorded as 
such, and volumes of product removed from the system which 
are not registered by measuring, devices or otherwise 
recorded, may be identified. 

Another advantage of the present invention is that 
discrete one-time unrecorded removals of product from the 
system may be distinguished, from continuous Losses 

consistent with leakage. 

Another advantage of the. present invention is that 
product leakage from all parts of the system, extending from 
the fill point to the point of discharge, may be identified 

and confirmed. 

Another advantage of the present invention is that 
secular, seasonal trends and repetitive special demands to 
provide short and long term estimates for demand of the 
product, and optimal reorder quantities and delivery 
schedules, may be identified and determined. 

Another advantage of the- present invention is that a 
better estimation of product volumes and volume changes may 
be obtained by calculating the differential of the volume 
function. 

Another advantage of the present invention is that a 
fuel access control system may be used to monitor dispensing 
-activity from and control access to the fluid storage 
system. 

A further advantage of the present invention is that 
all of the foregoing may be accomplished in a fully 
automated system that requires no human intervention, other 
than as an option available to the operator to enter 
quantities of material reportedly delivered for comparison 
with those computed. 
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Other features and advantages , af the invention will 
become apparent from the following detailed description, and 

from the claims. 

Brief Description of the Dra wings 
Fig. 1 is a schematic diagram of a facility 
including an underground tank storage system. 

Figs, 2, 3 and 4 are si portion of the Mathcad 
computer code used to perform the data, compression 
algorithm. 

Figs. 5, 6 and 7 are a block: diagram of .the streps 
performed during routine operation of the algorithm of the 
present invention . 

Fig. 8 is a block diagram of the steps performed 
during the data deletion operation of the algorithm of the 
15 present invention . 

Fig. 9 is .a block diagram of. the steps performed 
during the delivery calculation operation of the algorithm 
of the present invention. 

Fig. 10 is a schematic diagram of a data acquisition 
20 and transmission network that may be used in conjunction 
with the present invention. 

Fig. 11 is a schematic diagram of a facility 
including an above-ground tank storage system. 

Fig. 12 is a schematic diagram of a facility 
25 including a partially above -ground , tank storage system. 

Fig. 13 is a schematic diagram of an underground 
storage tank facility including a fuel access control unit. 
Description of the Pre ferred Embodiments 
The method and apparatus described herein applies to 
UST's, AST's or any type of storage tank. The product 
stored in the tank may be any fluid, including dry particles 
that flow in the manner of a fluid. 
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Fig. 1 shows- a UST facility- 1.0, illustrated as an 
automobile service station. Facility 10 includes a series 
of UST's 12, 14, 16 which may store - Che same or different 
types of liquid fuel product IS. Volumetric tank gauges 20, 
22, 24 in each tank measure the height of product 18 in the 
tank. Submersible pumps 26, 28, 30 in each tank pump 
product 18 to one of dispensing pumps 32, 34 through piping 
lines 36, 38, 40. Alternately, facility 10 may be an. AST 
facility with above-ground tank 1000, as. shown in.. Fig. 11, 
or a facility with a partially above-ground tank 1010, as 

shown in Fig. 12. 

Tank gauges 20, 22, 24 are mounted in tanks 12, 14, 
16. The tank gauges may consist of or be based on 
rnagnetostrictive tank probes or other sensing technologies. 
In the case of rnagnetostrictive technology, two floats 42, 
44 surround each probe, e.g., gauge 20 in tank. 12. One 
float 42 floats on the upper surface of product 18 in tank 
12, and the other float 44 floats on the interface of 
product 18 with any water or other foreign material 
collected at the bottom of tank 12- Tank gauge 20 
determines the distance between floats 42, 44 to obtain the 
height of product 18 in tank 12. Tank gauge 20 also 
contains temperature sensors 46, 48 , 50 spaced along its 
length to monitor the temperature of product 18 at various 

depth levels. 

Each of the dispensing pumps 32, 34 consists of a 
totalizer or flow meter 52, 54 disposed in a housing 56, 5 8 
to measure the volume of product 18 dispensed through hoses 
60, 62 and nozzles 64 y 66. To operate dispensing pump 32, 
nozzle 64 is removed front housing 56, which actuates 
dispensing pump 32 and causes product 18 to flow through 
hose 60 due to the pumping action of submersible pumps 26, 
28, 30 - a value , stored in totalizer • 52 is incremented as 



fuel is dispensed- through, hose 60 - Upon completion of the 
transaction, nozzle 64 is replaced in housing 56, thereby 
turning off dispensing pump 32 and discontinuing the action, 
of submersible pumps 26, 28, 30 and totalizer 52. 

Transactions are recorded electronically by software 
in a sales recording device 71 connected to totalizers 52, 
54 of dispensing pumps 32, 34. Totalizers 52 , 54 in 
dispensing pumps 32, 34 are connected to salsa, recording 
device 71 by means of communications and power, supply wires 
78, 80. 

Sales recording device 71 contains software capable 
of emulating the functions of a point of sale (POS) terminal 
associated with fuel sales made at facility 10. POS 
emulation software in sales recording device 71 functions on 
the basis of read only commands to eliminate the possibility 
of conflict with control commands from a POS terminal 
employed by facility 10. Alternative data acquisition 
systems can result in destruction of credit card. sales 
records, inadvertently shutting down the entire system, 
and/or causing electrical interference in the pump links. 

Tank gauges 20, 22, 24 are connected to a tank 
monitor 82 by means of communications, and power supply wires 
84, 86, 88 or communicate data through radio frequency 
transmission. Tank monitor 82 converts raw data obtained 
from tank gauges 20, 22, 24 into a form usable by a 
computer. 

A computer 70 contains a processor 72 capable of 
running various computer software applications and a memory 
74. Tank monitor 82 and sales recording device 71 are 
electrically connected to computer 70 to relay totalizer 
values, product height and temperature data to computer 70. 
Software executable by processor 72 of computer 70 is 
capable of querying- tank, monitor 82 .and sales recording 

- 12 - 



device 71. to obtain .measurement rdata at selected • time ■ 
intervals. The data is continuously evaluated as it is 
collected and is stored in memory 74 of computer 70 for 
later retrieval and detailed. ..analysis. r alternatively, 
computer 70 may communicate with a host processor 90 at a 
remote location. The continuous evaluations or detailed 
analysis may then be performed by host processor 90, which 
may be faster or more efficient than computer. 7Q. 

As an example, computer 70 may be a, personal 
computer or any other proprietary microprocessor-based unit. 
Computer 70 may capture data automatically through direct - 
connect serial interfaces with tank monitor 82 and sales 
recording device 71, or by manual operator keypad entry. 
Computer 70 communicates: with equipment at facility 10 
through four programmable serial communication ports, such 
as RS-232 communication ports. 

Computer 70 may, e.g., store tank dimensions- and 
product characteristics., and concurrent time and. date data 
along with the measurement data. Computer 70 may be used- to 
produce error and analysis reports as calculated by the 
software. It may. also have alarm event -initiated 
capabilities, such as when a leak is detected in any of the 
tanks. Such a computer system can accommodate facility and 
customer specific requirements while maintaining complete 
compatibility with other system components. 

The SIR method involves • reconciling volume data 
obtained from tank monitor 82 and volume data obtained from 
sales records. Sales transactions may be detected in a 
number of ways, including an electronic signal emitted from 
totalizers 52, 54, by .voltage sensing of control relays on 
pump dispensers 32, 34, or by observation of product removal 
using tank gauges 20, 22, 24. 
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It is essential that the- measurement used, to obtain 
these two types of data are made simultaneously. The SIR 
method of the present invention collects and analyzes 
observations of sales volume* and. can*, volumes which are 
derived simultaneously. Failure to collect both types of 
data simultaneously would bias estimates derived from 
separate volume measurements. ».„...„*<■ 
The SIR method properly accounts for. Che effects of . 
temperature, pressure and specific- gravity In 
product from two or more tan** may be blended, such as to 
achieve varying petroleum octane levels at pump d "^ ers 
32 34 When different fluid products are blended, the 
tanks are treated as one unit, and an additional parameter 
is introduced to determine the actual blend percentages. 

Data concerning the physical characteristics of the 
tank configurations and the accuracy of the ^-"W- 
and metering devices is collected during installation ^and a 
set-up phase of operation of facility 10 to create a . basis 
for subsequent statistical analysis. Information is then 
continuously collected so that the statistical °* 
SIR can be performed by computer 70 or host processor 90. 

Several procedures are used either singly or in 
combination to obtain the volume observ ations . »^ re 
the system configuration provides for determining whether 
hoses and dispensers associated with a given tank are 
active, the system is queried on a minute-by-minute basis 
or on he basis of another predetermined time 
determine the status of the dispensers, when 
dispensers are idle, the values from totalizers 52. 54. the 
taxi volumes (i.e. product heights in the tanks, and 
temperatures are recorded d ^ 

Second, submersible pumps 26, 28, 
determine on/o£f status, when it is determined that the 
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pumps are turned off, the values from totalizers 5 2, 54 are . 
read, and tank volumes and temperatures are recorded. 

Third, software algorithms used by computer 70 
detect and measure leads and/or lags between the recording 
of sales events and corresponding gauge and meter readings. 
When leads or lags are encountered and constitute a physical 
characteristic of the data measurement and recording system, 
constrained optimization, rather than unconstrained. .. 
optimization, may be used to determine parameter; estiraatres- 
Lagrange multipliers are one example of .such a constrained 

optimization method. 

The method of the present invention is capable of 
providing dynamic monitoring of system performance. For 
example, the leak detection function is carried out 
continuously while normal operations, e.g., removals and 
deliveries, are taking place. To detect leaks dynamically, 
the software is programmed to detect when sales or delxvery 
events occur and to calculate the volumes of product removed, 
or added as a result of such activities. Thus, dynamic 
testing does not require that the system be dormant and 
addresses the entire system from the point, of filling to the 

point of dispensing . 

The SIR method of the present invention also 
distinguishes between one-time removals and continuous 
losses consistent with leakage. The integrity or leak-free 
status of the system is not assumed a priori. Instead, the 
individual and unique characteristic pattern induced by each 
form of error when viewed along the separate dimensions of 
time, product height and sales volume are used to identify 
and quantify, the errors.. The method may also be used to 
detect and quantify undocumented removals, e.g., theft or 
additions of product. 
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Further, the overall system is self diagnosing in 
that it determines from the data the maximum degrees of 
reliability and precision of which a particular °P* rat ^ 
configuration is capable at any given.time, as well as the 
degree of calibration accuracy. 

in particular, product height in the tanks and 
temperature are measured continuously at, e.g., one-minute 
intervals . Height and gross volumes .are. converged to net 
volumes at, e.g./60-F or 15-C, using the algorithms 
described below. Sales recorded by the totalis 52 54 
are extracted and stored in memory 74 at times coincident 
with readings from tank gauges 20, 22, 24. If the 
dispensing system is capable of transmitting a signal 
indicating whether or not any or all individual hoses are 
active, that information is also stored in memory 74 
coincident with taking gauge and meter readings. 

The method of the present invention is designed to 
achieve the maximum accuracy possible, within the limi^tions 
imposed by the inherent random and irreducible noise m the 
various measuring devices incorporated. It ut ^"* a 
multiple measurements over extended time periods to identify 
and quantify systematic and repeatahle effects in the 
instrumentation and thereby correct for such effects using 

. . _ j: 4- U^-v Ae*tri r*<=X2 Trie 
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the known physical characteristics of the devices, 
system makes no a priori assumptions as to the ««»««*■ "t 
the devices used to measure product volume in the tank, to 
measure volumes removed, or as to the accuracy of volumes 
reported to have been delivered into the system. 

The resulting volumetric calculations are 
independent of the physical characteristics o£ the tank 
configuration and the various measuring devices which may be 
incorporated in the system. The results do not rely on 
npuTeacexed externally by the operator or from diagnostics 
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internal to the . measuring devices used:. Instead, the output, 
produced by the software which analyzes the measured data 
depends only on the patterns induced in inventory data 
produced by the tank gauges and measuring devices and, in 
5 particular, the cumulative variances that result when the 
various input values are combined. 

Various error patterns which the measuring devices 
can induce and. the effects of temperature, tank, geometry, 
and orientation on cumulative variances- are derived frorrr 

10 empirical analysis of real -world inventory data. The 

system's software synthesizes the output measurements of the 
various devices based on known characteristics derived from 
the empirical data. Thus, the software is capable of 
identifying measurement errors caused by the measuring 

15 devices and simultaneously compensating for the effects of 

those errors . 

Gauges can be systematically inaccurate in two ways. 
The height of the product in the tank, can he incorrect, and 
the height to volume conversion algorithms- may not reflect 
20 accurately the true dimensions of the tank or its 

orientation in the ground. The latter may be the result of 
incorrect measurements or an inappropriate conversion 
algorithm. 

The presence of such systematic effects and their 

2 5 nature may be established by examining the pattern of 

inventory variances as a function of product height. Errors 
of this kind induce patterns which repeat themselves as the 
tank is filled and emptied. If the tank length is 
incorrect, a linear pattern is induced. If product height 

3 0 is in error, a curvilinear pattern results, reflecting the 

varying volumes in different cross sections of a cylindrical 
tank. Tilt along the length of the tank induces a 
sinusoidal pattern symmetrical about the mid-height of the 
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tank. Absent such, errors, the .pattern will be purely 
random, reflecting only the inherent noise of the measuring 
devices. The absence of randomness and the presence of a 
systematic pattern serves to identify the presence of. 
systematic error. The pattern of a departure from random 
and its extent determines the source and extent of the 
systematic effects and the means necessary to correct them. 

Dispensing errors, unlike volume mwniring errors, 
are independent of product height, but are • sensitive to tfre 
volume of product dispensed. The nature and. extent of 
dispensing errors can be established by examining inventory 
variances as a function of sales volume. As in the case of 
volume measurements, in the absence of systematic errors-, 
variances as a function of sales volume will be random. The 
form and extent of departures from randomness serve to 
determine the source and extent of the errors and provide 

for their removal . 

Leakage from the system creates a continuous 
downward trend in the cumulative variance when viewed as- a 
function of time. By contrast, one-time additions and 
removals of product cause significant upward or downward 
translations of the cumulative variance which remain 
permanently in the record and do not introduce a continuous 
trend. Leakage is distinguishable from tank gauging errors 
when viewed as a function of product height because the 
pattern does not repeat as the tank is filled and emptied. 
If product is leaking from the system, a series of parallel 
translations in the cumulative variance is generated, each 
shifted by the volume of product lost between deliveries . 

The accuracy of measurements taken from the various 
components of the system determines the accuracy achievable 
in any one individual observation. Since the leak rate is 
computed from a. series of successive observations, however, 
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the minimum detectable .. leak rate can - b* reduced, to any 
desired magnitude by increasing the number of successive 
observations recorded. Thus, the system can serve as a 
final verification for leakage indications obtained, by. other 

methods . ' . 

At the conclusion of an initial set up period of 
data collection including one or more delivery and sales 
cycles , the collected measurement- data is analyzed, by 
regression analysis. The initial set-up regression is used., 
to derive tank dimensions and orientation, individual meter- 
calibrations and secular trends. A confidence level value p 
is computed at the .01 level of significance to determine 
the minimum leak rate detectable by the system, and the 
residual variance is computed to provide the current noise 

level of the system. 

The regression is performed according to the 

following equation: 



where : 
sti(R,L,T) 



Volume in gallons derived from the ath 
gauge reading in inches in a cylindrical 
tank with or without hemispherical end 
caps with radius R, length L, and tilt 
over its length of T inches. 

a = Initial inventory in gallons, which is to 
be estimated. 

Sa . Sales volume recorded on the kth 
kj totalizer. 

a = Fraction of sales volume recorded on the 
Jcth totalizer actually removed from the 
tank, which is to be estimated. 

D = Volume of the jth delivery. 
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Et = Elapsed time since . initiation until the 
ith gauge reading is recorded. 

Ls = Constant gain or loss in product per unit 
of time . 

B = Volume of product added (e.g. delivery) or 
Bj removed during some discrete time interval 

prior to or during observation period j. 



(0 if j < i 
r jiif jii 



All of the parameters are estimated simultaneously 
using least square estimation procedures. The R and T 
parameters are derived numerically, but the other parameters 
are derived analytically. 

Further, all of the parameters, including the 
initial inventory, are estimated simultaneously . The 
initial volume must be estimated from all succeeding data 
even if the tank is initially empty, otherwise the initial 
gauge reading and its conversion to gallons is assigned a 
credibility not assumed for all succeeding readings. Also, 
in a great majority of applications, the initial inventory 
in an already existing and operating system is not 

2 0 accurately known. . , 

initial inventory estimation is vital m determining 
the geometry of the tank. When tank geometry, tank 
orientation, or tank product height measurement depart from 
the values obtained from nominal sources, all gauge and 

25 meter measurements are affected. It is practically 
impossible to detect the errors induced in the gauge 
rneasurements and correct for them unless the estimation of 
the initial inventory is made coincident with the estimation 
of the values of the other parameters. 
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The estimate of the parameters are based, an the^ 
totality of the data collected. This means, e.g., that the 
estimate of leak rate Ls is determined from a linear trend 
including all of the data collected, not merely at one end • 
of the reconciliation period. Likewise, estimates of tank 
dimensions and orientation are derived from their overall 
contribution to reduction in residual variance, as opposed 
to a sale by sale analysis of tank, segments. 

The volume st^R^T) is derived from the. product- 
height measurement by multiplying the constant area of tank 
segments of height h (in inches) by tank length L. The 
volume in gallons of product in a horizontal cylindrical 
tank of radius R is given by: 



Vo1 = 231 



L R 2 cos -i - ( R - h) (2Rh - h 2 ) 2 



in the case of a tilted tank, the area of the 
segments varies with position, along the length of the tilted 
tank, and the volume is determined by integrating over the 
length L. Such integration does not result in a closed form 
because the cross sections are not circular, and a numerical 
integration would severely limit the frequency of 
observations. Instead, in this application the tank is 
treated as lying horizontally and the product is considered 
tilted, to derive an equivalent volume. This integration 
yields the closed form: 

Vol = _£L [ (z ~ l)8in-i {22 _ z a,T - 1 <2z - z*)^ Ml 
231 J 

The integrand is evaluated between the normalized product 
heights in inches., hu/R and hl/R, at the lower and higher 
ends of the tilted tank, respectively . It is standard 
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industry practice to install tanks on- an .incline to dxvert 
water and sludge away from the submersible pumps. 

Tank tilt is identified from the pattern it induces 
in the record of cumulative variances as a. function of 
product height . It is compensated for by fitting the 
correct mathematical form for height to volume conversions 
in a tilted tank to the cumulative variance calculated by 
the method of least squares. This is done simultaneously 
with estimation of the initial inventory. 

Tank length L and radius R are established by- 
equating the first partial derivatives of the sum of squared 
cumulative variance with respect to length and radius and 
determining the values which minimize the sum of squared 
variances. Simultaneous estimation of initial inventory is 
also required when estimating tank length I* and radius 

Errors in measurement of the product height h in the 
tank are characterized by curvilinear patterns induced by 
height to volume conversions in the cumulative 
a cylindrical container when heights are tra^sposed^ward 
or downward. Such errors also are compensated for by 
minimizing the sum of squared cumulative variances with 
respect to increments or decrements to measured product 
height. This estimation also requires simultaneous 
estimation of the initial inventory of the tank. 

in general, the accuracy of the estimates of the 
tank dimensions, tank orientation and height - a -™ tS ^ 
confirmed by observing that the cumulative variances of each 
derived value as a function of nominal product height are 
random and display no systematic influence or effects 

Dispenser totalizer calibration is continuously 
monitored and evaluated by minimizing the sum of squared 
cumulative variances with respect to multiplicative 
constanta associated, with, individual reported cumulative 
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sales volumes from all pump .dispensers associated with a 
particular tank system. This eliminates the need for manual 
verification of meter calibration. 

In particular, gauge performance is- continuously 
5 monitored to identify gauge malfunctions or degradation in 
gauge performance. Monitoring of gauge performance is 
independent of diagnostics which are internal to the 
measuring device. Diagnoses of problems ar:e based only on. 
their impact on the cumulative inventory variances, which are 
10 continuously monitored by the software. 

If the gauge fails to record changes in product 
height when the dispensers register sales, an increase in 
cumulative variances approximately equal to sales volume is 
observed; this effect can be identified by the monitoring 
software and a warning of gauge malfunction generated to the 
operator. 

However, observation of the gauge registering 
product height change, but with, a time lag after. sales are 
recorded, may be a feature of normal gauge performance. 
Such normal gauge performance is identified by repeated 
positive increments in cumulative variances as sales are 
completed with subsequent return of the cumulative variance 
to normal bounds. When such gauge function is determined to 
be the normal operating characteristic of a particular 
25 system, constrained optimization with lagged variables is 
introduced into the software. Otherwise, the gauge's 
performance is reported as a malfunction. 

Finally, temperatures in the tank are monitored to 
detect changes that are excessive for the time intervals 
between observations . Erratic temperature readings are 
deleted, and may indicate gauge malfunction. 

The software computes actual, rather than nominal, 
delivered quantities, and requires no input by the system 
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operator. The operator, may choose to input into the system 
the nominal delivery quantity indicated by the delivery 
invoice, along with the temperature and coefficient of 
expansion of the product- at the point of pick-up. The • 
software will then compute overages or shortages between the. 
nominal and actual quantities delivered, as well as the 
overages or shortages caused by temperature- induced 
variations in the transport of the product, to the .facility, 
and in the subsequent mixing of the. delivered product w±th. 
that resident in the tank. 

Delivery is identified by the software when a 
positive cumulative variance is observed which exceeds the 
system noise level and is not succeeded by a return to 
normal variance bounds. Delivered quantities are computed 
by estimating the volume increases they induce in multiple, 
successive observations. The required number of successive 
observations is determined as that sufficient to generate a 
confidence width which is within a predetermined tolerance. 
The system of the present invention is capable of accounting 
for sales conducted during delivery and for. noise introduced 
by post delivery turbulence in the tank. 

One-time unaccounted for removals or additions to 
the tank are computed in the same manner. Deliveries are 
distinguished from' such events by computing the rate of 
input, which in the case of normal gravity delivery should 
exceed 100 gallons per minute. . Other modes of delivery, 
e.g. pipeline delivery into above ground tanks, are 
identified by incorporating their known delivery rates. 

Leakage from the system is identified by a 
continuous linear negative, trend in the data which exceeds 
the computed minimum detectable leak rate after all of the 
various error phenomena described above have been identified 
and compensated for.. This calculation deals with the 
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totality , of the data obtained, by constantly monitoring known 
removals and is not restricted to observations made , only 
when the system is dormant. It is also independent of any 
single data reconciliation calculation in that trends 
throughout all of the data are evaluated. 

All calculations concerning volumes are made on the 
basis of net volumes, according to the following 
definitions: 

Net Volume , „ . 

in Tank = Gauge Volume (1 - (t-SO)CETJ 



where 



t - Measured temperature in degrees Fahrenheit 
(if centigrade, the term in parentheses 
becomes (t-15) ) . 

CE = Coefficient of expansion. 



and 

Net Sales 



Volume = Metered Sale (1 - ( £x-±_t* - 60) CE) 



2 



where tl and t 2 are temperatures measured by the tank gauge 
at the beginning and ending of a sale transaction, 
respectively. Deliveries are computed in net gallons, but 
are converted to gross quantities if required, based on 
external information input by the system operator, as 

follows: 

Gross gallons on invoice at the 
originating terminal. 

Net gallons on invoice at the terminal. 
tT = Temperature at the terminal. 
CE. . = Coefficient of. expansion. 
The program also records: 



GT 
NT 
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tA Ambient temperature in "..the tank, prior to ..: 

delivery. 

tF = Temperature in the tank at the conclusion 
of delivery. 

5 The following value is computed: 

tS » Temperature of the product in the delivery 
vehicle at the facility at the beginning 
of delivery. 

m . tF + NVA (tF - tA) 
10 NVD 



where 



NVD = Actual net volume delivered, previously 
computed . 



NVA = Net volume in the storage tank at the 
15 start of delivery. 

NS = Net overage ( + ) (underage (-)) in delivery. 
NT-NVD 

GVD = Gross volume delivered* 

NVD (1 -h (tF - 60) CE) 

20 GVS = Gross volume in the transport vehicle at 

the facility prior to delivery. 

NVD (1 +■ (tS - 60) CE) 

GSM = Shrinkage due to mixing in the tank. 

GVS - GVD 

25 GVT = Actual gross volume in the transport 

vehicle at the facility. 

«• - NVD (1 + (tT - 60) CE) 
GST * Shrinkag-e during transit to the facility. 
GVT - GVS 
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GOS = Gross overage (-<-) (underage (-) ) adjusted 
for temperature effects. 

GT - GVD + GST +■ GSM 
Calculations of volumes actually delivered are based 
5 on multiple observations of the balance of measured tank 

volumes and cumulative sales. This method requires frequent 
simultaneous observations of sales and. in- tarik volumes (i.e.- 
product heights) and temperatures- 

The volume of product in a. tank is derived by. 
10 measuring the height of the product and using the geometry 
of the tank, which is assumed to be known, to compute the 
corresponding volume. In many instances, tank, dimensions 
vary substantially from assumed design dimensions. 
Regulatory specifications permit up to 10% variation in 
15 length and diameter of cylindrical tanks. 

Tank orientation can also cause complications in. the 
calculations. The volume corresponding to a measured heigh.tr 
varies substantially when the tank is tilted away from 
horizontal or rolled away from vertical. 

Further, tanks may also fail to conform to a known 
geometry either through faulty manufacture or installation, , 
or may suffer significant deformation during the course of 
operations. For example, many fiberglass tanks sag or bend 

along their length. 

25 m addition, installed tanks are typically 

inaccessible, and difficult to measure. Thus, it is 
necessary to confirm the accuracy of height to volume ■ 
conversions from generated inventory data and to identify 
and correct discrepancies where they exist. 

3 0 The foregoing problems are compounded when two or 

more tanks are manifolded together. Manifolded tanks are 
. joined together by piping systems and serve common 
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dispensers.. Tuus. sales quantities from manifolded «nks 
constitute withdrawals from all tanks in the manifolded 
system, but not necessarily in equal quantit.es. Product 
heights typically vary from tank to tank, but tank . 
geometries^ dimensions and orientation may also vary so that 
a procedure for correcting height to volume conversion 
errors for a single tank will not apply- 

The different factors which influence inventory data 
manifest themselves in distinct ways which, facilitate Chen, 
identification and correction. These factors are most 
easily identified by examination of their effects on 
cumulative departures of actual measured invent ory from a 
theoretical or book value when viewed across a variety of 
dimensions. In particular, one-time undocumented physical 
additions or removals of product, e.g. over or under, 
deliveries and pilferage, are evidenced by an addition or 
subtraction of a constant quantity from the cumulative 
variance at the time of occurrence and all _ 
observations. Continuous loss of product accumulating over 
time. e.g. leakage, is evidenced by a loss 

continuous loss of product varying proportionally with, sales 
value, such a line leak or meter miscalibration, may be _ 
determined by identifying a constant negative trend that 
cumulative only over periods where delivery lines are 
pressurized.^ ^ ^ ^ ^ ^ _ 

cyclically as the tank is successiveXy filled and emptied 
with no long term gain or loss of product, is the pattern 
associated with height; to volume conversion error^ The 
pattern. is cyclical because the error source "/ d »"= al 
each cycle as the tank is filled or <«*>» ed . " " reCraces 
distinguishable from the other patterns m that it retrace 
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die same path without the translation- which, would occur if 
physical loss or gain of product were taking place. 

This problem is most readily diagnosed by analyzing 
cumulative variance as. a function, of product height. If the - 
variances are random with no evidence of systematic, effects, 
height to volume conversions may be assumed to be correct. 
If not, the form of the induced pattern indicates the nature 
of. the conversion error. Thus, an. error la -tank length 
induces a linear pattern, an error in tank tilt induces a 
sinusoidal pattern, and a constant error in tank height- 
measurement induces an arc-like pattern. When other sources 
of loss or gain are present, the conversion error patterns 
remain, but are translated in each succeeding 
filling/emptying cycle to reflect the physical loss of 
product which has occurred during that cycle. Thus, 
confusion between conversion errors and other effects can be 
elitninated. 

Sales readings and product height measurements must- 
be made simultaneously. Since the number of observations in 
any one sales cycle is typically too few to generate a 
conversion table of sufficient detail to be of practical 
use, subsequent sales cycles and their corresponding 
deliveries must be incorporated. If, however, deliveries 
are unmetered and are used to approximate the volume (as is 
the standard industry practice) , significant inconsistencies 
are introduced. If an overage or shortage occurs during 
delivery, then all subsequent sales volumes correspond to 
tank cross sections which- have been shifted upward or 
downward from their predecessors. Averaging or statistical 
treatment cannot overcome this deficiency since there is no 
means of knowing without metering whether, by how much, and 
in what direction the data has been shifted. 



„f rtu nrejent invention may include 
The procedure o£ the present. __.<«. 
• •„„ tf heioht to volume conversion error is a 
:Z L 'he error is a problem then the system must 

s TaL orientation, height measurement — - 

— - d whether the — ^ ^ an on . 
by other gams and losses I ^ no 

site leak detection investigation is a^^, 
leakage is indicated, and one or all ofi ^ new 

0 rve-rra^^^ — - 

encountered, the e.ct ^TSS ^ 

„„,,,„ dispensed from each dispensei 
IS actually dispens polynomial using a. 

physical measurement. R * s ^ co convert 

variable of measured product heigh < ^ S ^ ^ 

r^^-rre-direrCrb™ IsLeTproduct height 
20 :r;Xt L heginning and ending o £ saies events 
WhiCh nr. W^^ dispensed entities are 
regressed polynomial ^.^J^ST^.. 

ro easured P^^^HH ^tn tnt Polynomial is 
25 subject to the constraint tna diameter , the result is 

evaluated at a height ^ " ^« deli very 
the total tank volume. OBserv*^ 

events are discarded. # 
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A f if tli order polynomial has .proven adequate in .most 
cases. Residual analysis maybe used to determine adequacy 
of the polynomial in the presence of severe tank 
distortions, and higher order polynomials may be introduced 
5 as necessary. The number of observations required is 
determined by estimating a confidence bound around the 
resulting polynomial with a width adequate for the desired 

resolution. Thus, 

ASalei = Actual, dispensed volume in period i. 

10 n = Product height upon, conclusion of 

ASale^. 

h = Product height prior to commencement 

1-1 of ASale t and after completion of 

ASale^. 

15 a Diameter of- tank:. 

Vol = Total . volume of tank. 

The converted volume for height his then given by: 

Vol (h) = a t h + a 2 h 2 + -.-+- a„h B 
The omission of a constant term in. the regression implies 

2 0 that 

Vol (h) « 0 when h « 0 
This ensures that the polynomial derived from the height 
differences is well defined. 

For manifolded systems, actual sales are regressed 
25 simultaneously on individual polynomials based on the 
various height differences in the several tanks which 
correspond to a particular sales volume, subject to the 
constraint that each polynomial evaluated at the 
corresponding .tank diameter yields the total, volume of that 

3 0 tank. 
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ASale t - a.xdk.u. - ^) +- aiz^t-u. - 

+ . . . + '^(h^u - h" tl ) 

+ .aullk.u - ha) +• aa'^a " h2 " ) 
+ . . . + a n2 (h n t . 12 - h B u ) +■ . • • 
• 5 ■ + a XB (h 1 . l . - h im ) + aj.diVx. - h 2 im ) 

+ . . . + anntd^t-lm ~ 

where : 

ASalei = Actual Sales volume in period i . 

v, - Height of product in tank j after 

h *-" " completion of Asale^ and prior to 

10 commencing Asale A . 
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1, 2, 



m 



H = Height of product in tank j after 

nij completion of ASale t - 

Number of tanlcs manifolded- 

15 m 

volume conversion for the m measured Heigntsr, h,, h 2 , . - . 
h„ in the total system is: 

voHhi,h 2 hj = J2 J2 a a hi 

where: 

Height of product measured in the ith 
n i tan k in the manifold. 

Delivery inaccuracies have no impact on this 
calculation since all observations made during deliver! 
are discarded. Height changes are related only to tne 
corresponding volumes dispensed. 

Prior determination of actual quantities dispensed, 
as opposed, to metered quantities, ensures that the only 



es 
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remaining source- of .error is random measurement error . . 
Regression is designed to accommodate random error of this 
kind and to facilitate inferences when errors are present. 

An alternative method of estimating volume of 
product based on product height in single or manifolded . 
tanks involves determining a volume function by integrating 
a differential of the volume function. The total 
differential of the volume function is estimated using one 
of several procedures, e.g., least squares estimation. For 
example, for a manifolded system of storage tanks, a-£. 
Sa; = it h 2 t h mi ) - V(h litl ,h 2 ^ h m ^) 

where Sa t is the measured volume change associated witir 
measured changes in product height during a dispensing event 
from the manifolded tanks, then 

Sa, = V(h llt h 21 4.x> " V^^.tL,^ 4^ 

where Vj ( hli h. 4 ) is the partial derivative of the 

volume function with respect to h,, the height of the fluxd 
in the j th tank. The least squares technique provides 
maximum likelihood estimates because measurement errors 
occurring in tank gauges 20, 22, 24 have been established to 

be normally distributed. 

A differential function for a volume function having 
any functional form may be estimated in this manner. For 
example, a high order, polynomial may be used and constrained 
to have a preset volume at a maximum height, zero volume for 
zero height in all tanks and/or zero value of the first 
derivative at maximum height and at zero height. 
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For example^ if - . 

hij - product height in_ tank f^^cojpletion 



the start of sale 



of sale j and prior to 

Sa k = volume dispensed in sale Jc, 
and the volume function is an r th order polynomial in the 
form 

m r 

£ X 



voi (i^,^ 5- ati " &J 



then 



2«1 j«l 



10 



where the linear term of the polynomial is omitted to 
provide a zero derivative at h=0 . Then, the following 
equation may be minimized, 



subject to 



r-l 



and 



2 a.^ hmaxi = VolmaXi 



j*2 



for 



i=»-l ,2,.. 



- 34 - 



^^^^^^^^^^^^^ 



WO 01/01 W6 . 



pcnvuGoo/ifiTg 



vhere hmax, is the .maxiroum product . height in tank ± , Voltnax* 
is the preset maximum volume in tank i, and in is the number 
of tanks in the manifolded system. 

The foregoing equation works well for m=l. For m>l , > 
5 a further constraint is required to ensure upward concavity 
of the individual volume functions near zero volume. This 
is accomplished by constraining the second partial 
" derivatives of the individual volume functions, to be 
positive at zero volume. . In the case of polynomial 
10 . functions and tanks with equal radii, this reduces to the 

constraint a x x =a 2 i = * • * ~ a «t i' ■ _ 

Alternatively, the volume function may take the form 

' where £ <h> is a function of the height h. The 

V(h) is 

.... e^^f(h) 

15 Numerical minimization may be used to estimate this 

derivative function. An advantage of a function in the form 
of V (h) is that it asymptotically approaches zero (0) near 
h=0 and one (1) near the maximum height. 

Determining the volume function by integrating an 

20 estimated derivative of the volume function has many 

advantages. For example, the data used to estimate the 
derivative consists of discrete measurements of dispensed 
volumes and corresponding product height changes, which 
avoids introducing ambiguities and errors due to inaccurate 

25 calculations of deliveries of the product. The data does 
not need to be sequential, and data for periods during 
deliveries and post delivery turbulence may be discarded. 



- 35 - 



■ 1 ' CTAJGOO/nrm 



10 



15 



20 



25 



Because the only error sources, are. in .the metering devices . 
(for which calibration may be determined as described 
herein) and random errors of height measurement (the 
magnitude of which may be determined as described herein) - 
the error resulting from the height to volume conversion may 
be contained within acceptable limits. Further, the volume 
function derivative may be estimated accurately because the 
system can collect a large number of data points, which may 
be stored in a compressed format as described below, and 
because the system avoids delivery calculation errors. As 
in the case of calculating the volume of product m a single 
tank, the sales, volume and tank height measurements must 
take place simultaneously, the calibration of individual, 
meters must be monitored and recorded, and a large volume of 
data must be collected and recorded. 

With respect to temperature, the temperature of 
product delivered into a tank system almost invariably 
differs from the temperature of the product already xn the 
tank. Its addition has the effect of expanding or ^ 
contracting the volume of the combined product. Thxs change 
" volume can create the appearance of incorrecf —ns 
of the height to volume conversion, appear as leakage where 
none exists, or it can mask the existence of actual leakage. 

It is therefore preferable, and frequently 
essential, that all volumes, sales, deliveries and product 
in storage be converted to a common ^™ r ? 
analysis. TYP-ally 60"F (WO is chosen as the standard. 
The conversion is accomplished as follows: 
Net Volume = Gross Volume (1 - (t-60)CB> 



30 where: 

t 



Measured product temperature in degrees 
Fahrenheit • 



= coefficient of expansion. 
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A3 ahove, all calculations are. in net gallons of: . product - . 

A complication to the calculation may occur if the 
tank gauges 20, 22, 24 used to measure product volume are 
designed for static or dormant mode tank testing. . Such tank 
gauges detect leakage when the tank is taken out of service, 
in this case, product volume changes due to temperature 
changes during the course of a test must be accounted for. 

Further, as shown in Fig. 1, temperature sensors: 4fi. 
48 50 are located at different heights in tank 12. If -the 
level of product falls below a given temperature sensor, the 
corresponding weighted temperature measurement is dropped 
from the average temperature calculation, and a temperature 
jump and corresponding volume change may be observed when 
the net volume is calculated using the new weighted average 
of temperatures. If uncorrected, such repeated jumps xn the 
data would preclude further analysis of the data, for leak 
detection or the generation of height to volume conversions. 

The system of the present invention may be used to 
overcome these temperature related problems. Using. the. 

following definition, 

NDB = The net cumulative variance in .the 
inventory data at observation BT. 

then, 

NDB N = a(l - (t 0 - 60)CE) 

-E Sa^l - (t/ - 60)CE) 
i-i 

-V„(l-(t, - 60) CE) 



where: 

Gross initial inventory. 
. Temperature . of init ial product volume . 

Temperature of product at observation i. 
Gross volume sold in period i. 
, Measured gross volume in. tank at period S. 
Coefficient of expansion- 
Absent random error or leakage, and assuming no deluxe*- 
of product, then 



a 

ti 

Sai 

V N 
CE 



10 



and 



NDB N 



a(l - (t 0 ~ 60) CE) 



IS 



20 



E Sa t (1 

i-I 



( ti - 60)CE) 



- V M (1 - (t M - S0>CE) 
Therefore, if . temperature jump to temperature t* occurs at 
an observation N + 1, then 



NDB N< 



25 



V M (1 - (t N * 60)CE) 

- Sa Htl (1 - (t* - 60) CE) 

- V N+X (1 - (t* " 60) CE) 
V N (1 - (t M " 60)CE) 

- Sa N+1 (1 - (t* - 60) CE) 

- (V. - Saw) d " (t * ' 60)CE) 
V N (t*-t„)CE 
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When, tirts final quantity HnB*. x is added to the 
volume where the transition occurs between temperature 
sensors, and all subsequent volumes, the effect of the 
transition is eliminated, and analysis proceeds as it would , 
where individual temperature readings are available. 

A large number of variables must be estimated by the 
software to implement the SIR system of the present 
invention. For example, as many as forty hoses and 
independent totalizers per tank system, as well as • 
deliveries numbering four or more per day must be . 
accommodated. Thus, a very large volume of data must be 
accumulated, encompassing a substantial spread of sales 
volumes from each totalizer for both the set-up analysis and 
subsequent routine monitoring. To accommodate this volume 
of data within current or conceivable future practical 
computer memory capabilities, the algorithm implemented by 
the software utilizes a matrix formulation which invokes the 
property of a sufficient statistic to reduce the. memory 

requirement . 

The calculations used to determine the various 
error, loss trend and delivery estimates have the form: 

B = (x^l-'xV 

MSE = fy-yFO T (v-xB) 

m+1 

S 2 = (x T x)" 1 MSE 

where : 



B 
x 

30 y 



Column vector of m parameters to be 
estimated. 

Matrix of parameter coefficients. 
Column vector of independent variables 
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MSE = 



Mean squaxed. error- 

Variance covariance matrix of parameter 
estimates . 



The values contained in vector y comprise can* gauge 

„•„„. The entries in matrix x are measured sales 
readings. The entries ,_ lues The parameters of 

volumes, time, and other constant values^ T P 
vector B which are to be evaluated include ^J^^ 
volume of the system and subsequent volume cbanga*. 

including delivery amounts. „ mrfe< i .very- 

For example, if observations are recorded every 

minute, as many as !440 rows in the x matrix 

vector may be recorded. It would clearly * »*££Zt 

accumulate and store data in that form over an 

period of time. Instead, data compression techniques are 

Cited so that only a manageable amount of data need be 

stored. *-vraf if -an n x nt 

The algorithm utilizes the where 
matrix A is partitioned into two submatrxces. Band _ 
B is an i x m matrix and C is a J x m matrix, such that 
j « n, then 

s A T A +• B T B 



C T C 

For' example, _ 
^1 2 

if C - I 2 
2 1 
1 2 



then A = 



and B 



1 
2 



A A = 



1 
2 



1 
2 



1 
1 



2 
2 



J. L 



2 
4 



4 
8 
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B T B = 



2 
1 



Thus, A T A + B T B - 




-7 



S 4 
4 5 



C T C 



r 

2 



2 
1 



~1 



1 
1 
2 
1 



2 
2 

l' 
2 



7 
8 



a 

13 



At the conclusion of each 24 hour or other period, 
only x T x and x T y are computed and stored. The matrix x has 
the form of a square n x n matrix. Further, the aggregates 
of observations for different periods are additive, since 
two square matrices having n x n dimensions may be added. 
Thus the total data storage requirement for each period is 
determined only by the square of the number, of parameters of 
interest . 

The system is able to accommodate virtually 
unlimited numbers of observations by this method of data 
compression. Without this capability, the system would not 
have the storage capacity to accurately and simultaneously 
estimate the numbers of parameters which are required to. 
perform a statistically significant calculation. This data 
compression method also allows for processing the data at 
the facility or for transmitting the data to a host computer 
for periodic analysis. Figs. 2, 3 and 4 show the Mathcad 
computer code used to perform the data compression 
algorithm. 

Furthermore, (x T x)' l x T y is a complete and sufficient, 
statistic for B. No statistically useful information is 
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remaining is the precision 

used. -> - including 

The software performs SIR analyse, x 

ipak detection, using- ti«= 
inventory estimation and leak dec 

equation in the following form: 



10 where : 



and 



15 



x 

y 

range * 
meters = 

^ range ~ 
T range ~ 



, v <meters> 

\ ^ range ' 



(Strange) " (Grange) 

! . (number of observations) 

! ' . . (number of dispensers) 



Column of 1 ' s . 

Cumulative time in minutes. 

Cumulative sales for- individual 
dispenser in gallons. - 



20 



CD range 
S t JCj-ange 



25 



30 



Cumulative deliveries. 
. Tank stick reading in gallons ; 
• - initial inventory, ttut matrix rc 

To estimate the inlt ^ TQ . estimat e loss 

reS u 1C for Che 

B . Loss trend. 

Individual meter error. 
B 2 ^necers „,-r=niPter estimates, 

9 is the vector contain^ loss . 

namely beginning inventory, row {n «2) . 

rate.. The loss: rata-. estimate xn th 

- 42 - 



S 2 is the variance cavariance • matrix of-the parameter 

„ c _ /q2 \ 1/2 is the standard deviation of 
estimates. Thus, S 22 (S 22 > is cue 

the loss rate estimate. Finally, the minimal detectable 
leak is defined as t Aa , where t.. is, the U-«> percentile of. 
the Student's t distribution. 

The software performs delivery calculations using 
the equation in the following form: 



x 



[lrange / v x range 
y =r (Stk range ) 

where : 

range = 1 . . - (number of records) 

and 

i = Column of 1' s 

J- range 

T r 



L range 
^range 

^ range 



Cumulative time in minutes. 
Cumulative sales in gallons. 

0 where T raage is less than, daliveiry time 

and 

1 where T range is greater than or. equal to 
delivery time. 

stIc = Tank stick reading in gallons . 

The values of B, MSB and S 2 are then calculated, 
producing the following result for the vector B: 
Bl - Estimated initial inventory. 
3 2 _ Loss trend. 
Bj « Meter error. 
B 4 = Estimated delivery amount. 
S 2 is the variance covariance matrix of the 

_. e _ t&i )U2 j_ s the standard deviation of 
estimates. Thus, S 44 - (s 44 ; is we 

B the delivery volume estimate. The delivery tolerance is 
B^S*. -here t. is the U-a) percentile of the Student's t 
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■ ta& reduced, to" any . 
distribution. Delivery "^"^ ^ Rations used 
desire* value by increasing the renter 

in the calculation. r^chod. of the present 

T ne SIR analysis used * cumula tive 
invention involves computing ano^ari^^ cotnpu ted 
variances. when the initial set P ^ £orwa rd an 

tr end and meter calibrations are used t op ^ ^ 

ejected cumulative variance. ^ 

C ne difference between gauge £hen compu ted 

inventory. Actua! cumulative variances co 
from all subsequent gauge and meter reading 
Che expected variance. ope „ ti 

Figs. 5, 6 ana 7 D „ ftuare to perform this 

procedure 100 '^J? system and «- most 

S analysis. Data from the set up ^ ^ iM _ 

recent analysis is entered into the pg . <B-B . 10M , 
Tne data entered includes the erIor „, 

C anK tilt, meter ^^"^^i-a- established 
calculated trends. At step 10 . 3 , are seC at 

0 as counters. Counterl. Counter I and ^ 

zer o. The measurement data (rom ^ dtsp enser 

entered at step 1... ^'f^the product temperature, 
totalizers, the product height ^ th • P ^ ^ 

The software computes the 9 rass cbe sal es as 
(S product, the most recent gross volume and^the ^ ^ 
measured by the individual a«P en 

110 by 

software further manipulates ^ ^ 3 _ computlng 

converting all gross volumes to net ^ 
obs ervation « ^ervation^ariance^ is 
30 cumulative variance. The. sign 

recorded at step 112. cumulative 
The program proceeds 114 , 120 , 124. 

fiance and the value of Counterl » steps 
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128 132, 136 and 140 : Depending on tire, cumulative variance, 
and' the value of Counter!, the program analyzes the 
collected data at step 118 if it is a final observation 
(step 116), deletes the collected data (steps 122 and 134) • 
performs the analysis for a delivery (step 126) (see below), 
or reads new data (steps 116, 122, 130, 134, 138 and 142) 
upon updating the value of Counterl and other computational 
variables (i.e. index, sign index and sign) . In .some cases , 
collected data is deleted (steps !22 and. 134). 

Upon computing the loss rate at step 144. the. 
program reads new data at step 146 if the loss rate is not 
greater than or equal to, e.g., .2 gallon per hour, 
otherwise it computes the trend of the data at step 148 
at step 150 it is determined that the trend is greater than 
.2 gallon per hour, a warning is issued at step 156 In 
either case, the software continues to read, and analyze the 
data at steps 152, 154, 158 and 160 until, the last 

observation- 

The operation of deleting data 170 is shown, isi 
detail in Fig. 8. After performing similar analyses at 
steps 172, 174, 17S. 180 and 184, using the indues and the 
values of the calculated standard deviations as in the 
routine operation procedure described above, the values of 
the counters are updated and new data is read at steps 178 
182, 186 and 188. Data is deleted in accordance with steps 

17a, 186 and 188. 

Finally, Fig. 9 shows the delivery calculation 190 
in detail. After determining that the cumulative variance 
is greater than a predetermined value (three "andard 
deviations) at step 192, the program determines whether the 

^, „ a „ -rnn aallons per minute (step 
variance is greater than, e.g., loo gaxion* p 

194) If so, the delivery is recorded and the amount 
delivered is determined at steps 202 through 220 . 
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If there is a delivery in. progress (step 20.2) . data 
is read until a negative observation to observation valance 
s observed <step 204, . The variance is monitored^ until the 
turbulence in the tan, subside. (step _ 
5 observations are read (step 20S>, and a 

IS minutes before the delivery until the end of the 
turbulence observations are deleted (step 210) . An ^ 
indicator variable is introduced . with the ««^«^ 
observation, from which region If^,^' " 

10 The confidence bound on the indicator is computed (step 

214, If the confidence bound is within a Predetermined ^ 
tolerance, the volume of the delivery is ^ 
confidence bounds (step 220, , otherwise, ^ditronal 
observations are added, and the confidence bound is 

15 recomputed (step 218). less 
If the variance between data measurements- im ^ 

„ software determines whether 

than 100 gallons per minute, the software 
. the gauges are inoperative and reports them as being 
Inoperative (step l,s> . or procee a, as in the rout^ ±- 
20 operation procedure according to step 200 (in wh, 

a negative variance, depending J^^L oM standard 
variance exceeds a predetermined, value (within 

glelT Served variances are within , three 
25 standa rd deviations ^ 

;r "ire":; ::=:- *~ — * — 

Tother predetermined «^^^"1^c-. 
programs are executed depending on the nature a 

30 of the departure. (M Qther ^ 

observations after the initial departure, the 
variance. re.^ to within the tolerance range, all data 
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from and including the initial .departure and prior to the 
initial observation are deleted. The time extent and number 
of observations involved is recorded and stored for, e.g., a 
daily gauge performance report; 

If all ten (or other predetermined) successive 
observations remain outside the tolerance bound and the 
cumulative variances are of the same sign, a new trend line 
is initiated at the point of initial departure. After ten- 
tor other predetermined) additional observations, • a third 
trend line is initiated. If the increment to the overall 
trend estimated from the most recent observations is not 
significant, the most recent data is consolidated with. the 
previous data and the process is repeated until such tune, 
if ever, that the current trend increment is significant. 

If the departure is positive, the system checks 
whether the product is being dispensed and whether the gauge 
height fails to decrease, reflecting- removal from -tank. If 
so, the tank gauge is reported to be inoperative. 

If the gauge height is increasing-, tnonxtoring- xs 
continued as above until the most recent trend line returns 
to its original slope. Minute to minute variances are 
monitored to detect turbulence until the gauge values again 
return to within tolerance. All observations which occurred 
in the fifteen minutes prior to the first positive departure 
until the end of post delivery turbulence are deleted. An 
indicator variable is introduced at the first observation 
after post delivery turbulence. The system collects thirty 
additional observations and performs the regression from the 
beginning of the period to determine the volume delivered. 
3 0 The volume delivered is then reported. 

If the departure is negative, the system proceeds as 
with delivery. If successive slope increments fail to show 
a return to the original slope, indicating continuing loss. 
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of product for -.a predetermined period, typically one hour 
and the slope exceeds .2 gallon per hour, the system reports 
a warning that there is a continuous loss of product. If 
the loss rate is less than .2 gallon per hour, but greater 
than the minimum detectable leak, the system continues to 
monitor and recalculate the parameters, to be included in a 
daily operational report. If the incremental trend ^™ 

shows a return to the original trend, the system proceeds, as 
with delivery, introduces an indicator variable, deletes 

data as necessary, and performs the regression to determine 

the volume of product removed. The system reports a one- 

time removal of product. 

Referring to Fig. 10. the invention incorporates a 

data acquisition and transmission network (DAT network) 300 

to completely automate the process of 

transferring and processing product inventory data f or use 
in product management. , delivery scheduling and envxronmental 
compliance practices. DAT network 300 includes 
processors 302. 304 at the facilities 306. 308 where the 
tanks are located, a customer host processor 310 and a 

ii -> hit network 300 links, multiple- 
central host processor 312. DAT networ* 

remote facilities 306, 308 to central host processor 312, 
which performs the SIR analysis. The link may be 
accomplished indirectly through customer host processor 310. 
which itself is connected to a plurality of 
facilities 306. 308. Each of these processor 
composed of independently operating software 
systems which form the basis of a wide area network linked 
by modems which transmit information electronically via the 
telephone network 314 using standard dial-up voice grade 
telephone lines. Examples of DAT networks are <*»J"^ 
and ECCOSIRA systems developed by Warren Rogers Associates. 
Inc., Hiddletown.,. Rhode Island. 
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DAT network 300 provides a-, uniform method, of 
integrated management for the widest possible variation of 
underground and above-ground fuel storage, movement and 
measurement systems. On-site processors 302, 304 are 
capable of obtaining information from any electronic or 
mechanical control system, enabling DAT network 300 to 
accommodate facility configurations that are unique to each 
facility while presenting the inf ormation captured at remote 
facilities 306, 308 to customer host processor 310 or 
central host processor 312 in a uniform format. 

On-line processors 302, 3 04 obtain and capture 
product inventory data through the use of proprietary 
interfaces with external systems in use at remote facility 
306, 308, such as tank gauges and sales recording devices. 
On-line processors 302, 304 transfer captured information 
daily, weekly or monthly through the public switched 
telephone network 314 to customer host processor 310 or. 
central host processor 312 for use in inventory management, 
delivery scheduling and/or environmental compliance. On- 
site processors 302, 304 may be, e.g., touch-tone telephones 
acting as sending units and Windows-based multi-line, voice, 
prompt /response PC's as the receiving units. On-site 
processors 302, 304 may be designed to meet the specific 
needs of facilities 306, 308 without requiring remote 
hardware at the facility in addition to that already 
present . 

in particular, each of on-site processors 302, 304 
typically is equipped with an alphanumeric keypad, a 
character display, a power supply, four programmable serial 
communication ports, an .internal auto -dial /auto -answer 
(AD/AA) modem and a local printer port (for connection to a 
printer) . The keypad and display allow for operator 
configuration and manual entry of sales, delivery and tank 
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Level data . Use - of . an AD/AA. ZAOQ baud- modeni allows multiple 
on-site processor 302, 304 to share an existing- voice grade 
telephone line by establishing communication windows to 
minimize attempted simultaneous use. Each of the 
5 programmable serial communication ports is independent, 

fully programmable and governed by options selected at the 
facility or off -site through modem access. Finally, on-site 
processors 302, 304 can prompt the facility operator- to 
enter missing or suspect entries when results are outside 

10 the expected range. 

The use of customer host processor 310, which is 
capable of receiving, storing and processing information 
from multiple on-site processors 302, 304, enables the 
management of a remote tank population from a single point 
of contact . A database of information created by customer 
host processor 310 is the basis for all higher level product 
management functions performed by DAT network. 3 00. The 
database is also the basis for the environmental . compliance 
analysis performed by central host processor .312- 

The use of central host processor 312, which is 
capable of receiving, storing and processing the information 
in the database created by customer host processor 310 for 
product management enables DAT network 300 to achieve 
maximum results by utilizing the database for environmental 
compliance without additional remote facility information or 
communication. Central host processor 312 is capable of 
transmitting a resulting database of the environmental 
analysis back to customer host processor 310 for printing 
and other customer record-keeping requirements. 

The processor elements of DAT network. 300 may 
exhibit other useful operational characteristics. .To 
prevent unauthorized access to DAT network 300, a security 
access code for dial-up. data transfer functions is required.. 
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Under secured access, the baud, rate, parity, stop. bit', 
parameters and communication protocol are determined at any 
of on-site processors 302, 304, customer host processor 310 
or central host processor 312. 

Another function of DAT network 300 is to monitor 
tank contents generally. DAT network 300 can be programmed 
to activate, e.g., an audible and visual alarm if the water 
level in the tank is too high (e.g., greater, than. 2 inches) , 
if the product level in the tank is too high (e.g-. , more 
than 90% of tank capacity) or too low (e.g., less than 10% 
of tank capacity, more product must be reordered, or less 
than two days supply) , and if a theft occurs (product level 
changes during quiet periods) . 

The system may be used to obtain valuable 
information other than inventory regulation and leak 
detection. For example, the system may incorporate time 
series analysis routines, including Box Jenkins, moving 
average and exponential smoothing, to derive estimates of 
demand for the product which also incorporate temporal and , 
seasonal trends and special events. 

The demand analysis may also be combined with 
additional inputs of holding costs, reorder costs, 
transportation costs and penalty costs for running out of 
stock. The system can include optimal inventory algorithms 
to determine optimal order quantities, reorder points and 
optimal delivery truck routing.. Further, the system may 
incorporate multiechelon, optimal inventory procedures to 
accommodate combined wholesale and retail operations, such 
as with calculus -based optimization and linear, nonlinear 
and dynamic programming . 

As shown in Fig. 13, a DAT network may include a 
fuel access control unit or system 510 at a storage tank 
facility 5Q0 such. as. a- UST automobile, fueling facility. 
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Fuel access control uait.Sia is- a. dispensing system actuated 
by the use of a device coded with information, e.g.. a card 
520 with a coded magnetic stripe 522, e.g. .an optxcal 
punched card, an electrically erasable programmable read^ 
only memory ( EEPROM ) key, a radio frequency identification 
(RFID) tag, a magnetic resonance coupler, a bar code, or 
other type of coded medium which contains identification 
information pertaining to the. user. Euel.access control 
unit. 510 may include apparatus for- a user to input 
information, e.g. , a card reader 512, a display 514, and a 
keypad 516, a control system 562 for turning a fueling 
dispenser 560 on and off, and a processor 564 or 
computing platform for controlling and monitoring the user s 
fueling process. Manufacturers of fuel access control 
system which rely upon optical reading devices or magnetic 
stripe card reading devices to identify the user include 
FillRite, Fuel Master, Gasboy, PetroVend and Trak 

Engineering. monitor the 

Fuel access control unit 5X0 is useu 
activity of fueling dispenser 560. Fueling dispenser 560 
includes a hose 566 for dispensing fuel from a tank SIS, a 
totalis 568 and a meter 569 for measuring the volume of 
fuel dispensed by hose 566. Fuel access control unxt 510 
may communicate with an on-site processor 530 located .ns.de 
facility building 540 over a local area network (L*K> . The 
ommunications between fuel access control unit S10 and on- 
site processor 530 may be over RS-232/KS-485/RS-48S 

(MulriDrop) cabling 542. 

Fuel access control unit 510 provides a system o. 
controlling access to fueling facility 500 by determining 
"^entity of each user of the f a=ility and screening each 
user based on his or her authority to purchase fuel, 
ratification of the- user is made by requiring the user to 
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present a valid magnetic card. (e.g. , card. SZCU , am . optical 
punched card, an EE PROM key, an RFXD tag, a magnetic 
resonance coupler, a bar code, or other type of. coded 
medium, which contains identification information pertaining 
to the user. Such fuel access control systems are referred 
to as island control unit's or cardlock system. Further, the 
user may be required to present additional identifying data 
by other means such as buttons, key switches, or by entering 
information on keypad 516. Once the identification data is 
collected, fuel access control unit 510 determines the 
user's fueling privileges, and based on this information 
will either allow or deny fueling by the user. If fueling 
is allowed, fuel access control unit 510 will enable 
dispensing pump 560 for that particular user and monitor the 
fueling process. At the completion of the fueling process, 
fuel access control unit 510 will record the amount of the 
fueling transaction in a memory 567 and retain the recorded 
information for further accounting of the transaction. 

Fuel access control unit 510 may be used to perform 
a variety of functions, including the following: 

1. Identifying the user by reading a card or other 
coded medium and collecting the user's 
identification information such as a driver license 
number or other personal data. 

2. Collecting other pertinent data for analysis, such 
as an identification of the user's vehicle, the 
vehicle's odometer reading, a trip number, the 
trailer hub counter, the engine hour reading and/or 
a refrigerator unit, hour reading. 

3. Making authorization decisions, to determine whether 
the identified user is permitted to obtain fuel. 

4. . Enabling fueling by. enabling the .proper dispensing 

pump for the user. 
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. . . ^ lirrcr t±ie maximum, amoimt 

5. Monitoring f ueling by controlling- . w» 

dispensed. 

6 Turning off the dispensing system if no fuel is 

dispensed for a predetermined- period o£ time. 
5 7 . Recording the fueling transaction by storing the 

final amount of fuel dispensed. 

Reporting the fueling transaction to a processing 
location for inventory analysis or other analysis. 
There are two types of authorization- procures fiqr. 
.0 determining whether an identified user iS ^^^^ 
obtain fuel from a fueling facility. Fuel access control 
unit 510 .ay use either or both of these ^^^^ 
procedures. For the first method of ^^^'JX 
Authorization, fuel access control unit 510 collects the 
autnorizaci™. collected information ta 

15 user's information and forwards the collected inr 
an outside agent to make a final decision, aa t= 
not the identified user is permitted to obtain fuel from 
fueling facility 500. The. outside agent may . 
approval, along with fueling parameters - 
20 amount) . or a denial. Fuel access control unit 510jill . 
then inform the user whether or not fuel may be obtained. 
The outside agent may be. connected to fuelaccess contro, . 
unit 510 via a dial-up telephone line, a LBN or a direct 

communication link. , harT1fl1 

2S For the second method of authorization internal 

authorization, fuel access control unit 510 collects the 
user's information and compares the collected -formation to 
a data table stored locally to make the final decision a to 
„he"er or not to allow fueling. The locally stored table 
30 lay return an approval, along with fueling ^ 

= denial Fuel access control unit 
a maximum amount) , or a denial . ru 

510 will then inform the user whether or not fuel may be 
obtain The locally stored table may be housed directly 
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in fuel access. control.unif.5LQ,. in a controL device,a C the,, 
fueling facility such as orAite- processor 530 or carried on 
the access medium (e.g., card 520) used to request fueling 
authorization. The locally ..stored, table may also be ,- , 
imbedded directly in fuel access control unit 510 or 
accessed via a LAN inside the fueling facility's building 

540. 

Fuel access control .unit. 510 functions as an. 
additional point of sale: (POST device; similar to- saXes 
recording device 71 (Fig. 1). Fuel access, .control unxf 5XQ, 
responds to requests for hose status and totalizer and meter 
values in the same manner as a POS device. Fuel access 
control unit 510 also monitors each hose 566 and tracks 
status changes in the hose, including indications that the 
hose is idle, that a request for access is in. process, that., 
the use of hose 566 has been- authorized, . that .the base, has , 
been taken off its hook, that dispensing-, pump 560 xsv 
dispensing fuel with hose 5.66, removed from its. hook. .. and. ; . 
that the dispensing pump has been turned-, of f : and the: hose- xs 
idle again. 

Each detailed transaction: that- is completed by fuel, 
access control unit 510 may be retrieved by on-site 
processor 530 from memory '567 upon completion- of the 
transaction. The transaction information may be stored in 
processor 530 for further, analysis,. Further, based on the 
stored, detailed transaction information, a detailed site 
dispensing audit can be performed. Such a site dispensing 
audit would determine whether the volume claimed to be 
dispensed by fuel access control unit 510 actually ^ 
represents the. volume change in the UST or AST during the - 
same period as calculated by on-site processor 530 . 

in conventional cardlock applications as well as 
other transaction authorization, procedures , the processing 
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fuel access control' system ^j^^ic tan* 

site processor 530 that receiv. =s data fro. ^ 

gaug e 580 may also "rac"^ By contrast, an 

conventional "^uT act s c^ un^ ,10 « 
enhanced, integrated fuel aCC " S j\ e rt , tose ss(r ., snch. 

includ e ^ accurate ana V*^^, ray ^ ^ 

and enabied by authorization control ■ ■ ^ 

^ t-o be dispensing fuel. From cue e 
appearing not to be ai p h&ve sinip _ ly 

of fuel access control unit 510, the u ^ 
ranged his mind about purchas^el^ ^ ta . 

^■p nn-qite processor ,5iOi <* v"**--. 
perspective of on-site p br analyzing the. 

Lde about the dispensing . pump's M £ activity . 

tan, activity and °-*"^» ^" ^ her . . althou^ other. 
of totalizer 5S and ^ ^ cha same 

r^odHn sit! "/ocesso/sso may tract s hose a«i~ty 
independently - ^ ch e scope of -the 
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A method of monitoring a fluid storage and 
dispensing system including a dispensing apparatus, 
the method comprising: 

dispensing a quantity of fluid from the system 
using the dispensing apparatus based on an 
authorisation; 

collecting a plurality of measurement data from 
the dispensing apparatus in a format readable by the 
computer; 

collecting a plurality of dispensing apparatus 
activity data including a state of the dispensing 
apparatus in a format readable by a computers- 
storing the plurality of measurement data and 
the plurality of dispensing apparatus activity data 
in a memory; 

determining the mathematical relationship 
between a volume of fluid in the storage tank and 
height of fluid in the storage tank; and 

statistically analysing the stored plurality of 
measurement data and the mathematical relationship 
to calculate a volume of fluid based on the 
plurality of measurement data collected from the 
dispensing apparatus. 
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The method of claim 1 wherein the authorisation is 
based on information relating to a user requesting 
that the quantity of fluid be dispensed from the 
system . 



The method of claim 1 wherein the authorisation as 
provided by a fuel access control system. 
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The 



method of claim 



further 



comprising 
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transmitting the plurality of measurement data to a 
processor to perform the statistically analysing 
step. 

5 5. The method of claim 1 wherein the dispensing 
apparatus includes a totaliser and a metre. 

6 The method of claim 1 wherein the plurality of 
measurement data is stored in compressed matrix 

10 form. 

7 The method of claim 6 wherein the compressed matrix 
is a product of the data matrix and a transpose of 
the data matrix. 

15 8 A method of monitoring a fluid storage and 
dispensing system including a plurality of 
dispensing apparatuses, the method comprising: 

dispensing a quantity of fluid from the system 
using the plurality of dispensing apparatuses based 
on an authorisation; 

simultaneously collecting measurement data and 
dispensing apparatus activity data from the 
plurality of dispensing apparatuses in a form a 
readable by a computer, wherein the dispensing 
apparatus activity data includes the state of the 
dispensing apparatus; 

repeating the collecting step to obtain a 
plurality of the measurement data; 

storing the plurality of the measurement data 

in a memory; and 

statistically analysing the stored plurality of 
measurement data to calculate a volume of fluid 
based on the plurality of measurement data collected 
35 from the dispensing apparatus. 
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A method of claim 8 wherein the authorisation is 
based on information relating to a user requesting 
that the quantity of fluid be dispensed from the 
system. 

A method of claim 8 wherein the authorisation is 
provided by a fuel access control system. 

The method of claim 8 further comprising 
transmitting the plurality of measurement data to a 
processor to perform the statistically analysing 
step. 

The method of claim 8 wherein the plurality of 
measurement data is stored in compressed matrix 
format . 

The method of claim 12 wherein the compressed matrix 
is a product of a data matrix and a transpose of the 
data matrix. 

A fluid storage and dispensing system, comprising: 
a tank for storing fluid; 

a dispenser for dispensing fluid from the tank; 

a fuel access control system for enabling the 
dispenser based on an authorisation; 

measurement apparatus for collecting 

measurement data corresponding to a quantity of 
fluid dispensed from the tank by the dispenser and 
for collecting dispensing apparatus activity data 
including a state of the dispensing apparatus; and 

a processor for statistically analysing a 
measurement data to calculate a volume of fluid in 
the system based on a mathematical relationship 
between a volume of fluid in the storage tank and a 
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ight of the fluid in the storage tank. 



15. A system of claim 14 further comprising a memory for 
storing the collected measurements data. 

5 

16. The method of claim 14 wherein the fuel access 
system is activated by a coded medium which contains 
identification information pertaining to the user 
and which constitutes a request to dispense fluid 

10 from the system. 

17. A system of claim 14 wherein the authorisation is 
based on information about a user requesting to 
dispense fluid from the system. 

15 

18. The method of claim 14 wherein the measurement 
apparatus comprises a totaliser and a metre. 



19. The system of claim 14 wherein the processor stores 
a measurement data in a compressed matrix format. 

20. The system of claim 19 wherein the compressed data 
matrix is a product of a data matrix and a 
transposer for data matrix. 

25 

21. The method of claim 1 wherein the fluid storage and 
dispensing system is a manifolded system. 
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22. The system of claim 14 wherein the fluid storage and 
dispensing system is a manifolded system. 
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